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Orthogonal polynomials on the unit circle, g-Gamma 
weights, and discrete Painleve equations. 



O 

psj Philippe Biane 



Abstract. We consider orthogonal polynomials on the unit circle with re- 
spect to a weight which is a quotient of g-gamma functions. We show that 
the Verblunsky coefficients of those polynomials satisfy discrete Painleve equa- 
tions, in a Lax form, which correspond to an A^^' surface in Sakai's classifica- 
tion. 



Jo 



cos{2TTxy)f{y)dy. 



< 

1- Introduction 

, ^ , Motivated by number theoretic considerations, J.F. Burnol [4], [5], has raised 

the problem of realizing the Fourier transform as a scattering. More precisely, 
consider the Fourier transform restricted to even functions on R which is given, 
^ after appropriate normalization, by 

On 

' One can rewrite this as an isometry of the space L'^{]0,+oo[,dx/x), by putting 

f{x) = g{x)/y/x, to obtain 

O Tg{x)^2l y/xycos{2Trxy)g{y)dy/y 

> . 

• ^ The interaction between the additive Fourier transform and the Mellin transform 

rS is captured by the formula, where the Fourier transform acts on tempered distri- 

^ butions, 

(1.1) :Fx -vr p(i/4^^^/2)'^ ''''^^^ 

which yields the gamma factor coming from the place at infinity in the functional 
equation for Riemann zeta function. Gamma factors for other L-functions are 
obtained by similar considerations involving the radial Fourier transform in R". 

Interpreting x**^ as an incoming wave, and a;"**^ as the corresponding outgoing 
wave (in logarithmic coordinates) on the space ]0,-|-cx)[, it is a natural question to 



ask what is the scattering potential which realizes the map (1.1) (we shall recall 
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the basics of scattering theory in the next section) . This was answered by Burnol 
[6] who showed that, more generaUy, the function 

r(a-zfc/2)r(6 + zfc/2) 
^ ■ ' r{a + ik/2)T{b-ik/2) 

{a,b real, a > b), is the scattering phase shift for a Dirac equation on ]0, +oo[, 
where the potential is given in terms of Painleve transcendents. In this paper we 
consider a discrete analogue of this problem. Discrete scattering problems have been 
introduced by Case and Kac [7j and studied further by Case and Geronimo |8]. In 
this last paper it is shown that, given a function /, analytic inside the unit disk, 
continuous on the unit circle, then /(e'^)//(e'^) is the phase shift for a discrete 
scattering problem, and solving the inverse problem for this scattering problem 
is equivalent to finding the recursion relations for the orthogonal polynomials on 
the unit circle with respect to the measure d9/\f{e^^)\'^. Recall that, if (f>n is the 
sequence of these polynomials, and 

then the recursion relation, due to Szego, can be written, with 



as 

(1.3) $„+i(z) = B„(z)$„(z) 
where 

(1.4) B„(z) 
The a„ are the Verblunsky coefficients of the sequence (j>n, see e.g. |21) . |13j . 



z an+i 



or [19j . |20j . Equation (1.3) is a discrete analogue of a Dirac system leading to 
scattering. Thus a natural problem, in view of Burnol's question, is to find the 
recursion coefficients for a weight w, expressed in terms of g-Gamma functions, of 
the form 



w{z) 



{bz,q)oo{bz i,q)c 



with, as usual. 





Indeed u;(e'^) = l/|/(e*^)|2 with /(e^^) = (6e*^ g)oo/(ae'^ gU, and taking a suit- 



able limit of /(e*^)//(e*«) for g ^ 1 yields the phase shift ( |l.2 |. 
The weight w satisfies a g-difference equation 

where V and W are two polynomials of degree 2. It turns out that a method, which 
goes back to Laguerre, further used by Freud 1_QJ, allows, in this situation, to obtain 
nonlinear relations between these recursion coefficients. This kind of question has 
seen renewed interest recently, see e.g. |17j . |22) . We shall prove that 

^n{qz) = :^A,,{z)^,,{z) 
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where An is a matrix with polynomial entries, of degree 2, which can be computed 
from the coefficients of the orthogonal polynomials. Compatibility between this 



equation and (1.3) implies the relation 

An+liz)Bn{z) = Bn{qz)An{z). 

This last equation can be interpreted as an isomonodromy deformation of the q 
difference operator associated with the matrix A„. It is similar, although it is not 
equivalent to, the Lax form of the discrete Painleve VI equation studied by Jimbo 
and Sakai in |15j . We shall study it in details and prove that the equation that we 
obtain in this way is a Lax form for a discrete Painleve equation, corresponding to 
the ^3^^ surface, according to the classification of Sakai |18j . As we shall see, this 
equation looks like the one of |15| . but is obtained from translation in a different 
direction in the afRne Coxeter group oi^K 

This paper is organized as follows. We start by recalling some basic material 
on one-dimensional scattering theory, following [9] , for the continuous case, and [8] 
for the discrete case. We also recall some of the theory of orthogonal polynomials 
on the unit circle, following closely Laguerre's method is used in section 3 to 
derive nonlinear equations satisfied by the recursion coefficients of these orthogonal 
polynomials. Then, in section 4, we give an explicit formula for these nonlinear 
relations as a birational map in two variables, and show that they give rise to a 
regular isomorphism between two surfaces obtained by blowing up x at eight 
points. In section 5 we recall Sakai's theory, and use it to identify our nonlinear 
relation as a discrete Painleve equation. We determine its type in Sakai's classifi- 
cation, and find its expression as a certain translation in the afEne Weyl group of 
type Finally in the last section we derive a differential system by taking a 

continuous limit of the discrete equations. 

The results of this paper are based on numerous lengthy computations. In order 
not to bore the reader we have omitted the details of most of these computations, 
but given enough details so that they can be reproduced. 

I would like to thank Jean-Frangois Burnol for keeping me informed about his 
work which motivated me to undertake the study of these equations. I thank also 
Alexei Borodin and Philip Boalch for their useful comments at an early stage of 
this work. 

2. One-dimensional scattering, continuous and discrete 

2.1. Scattering for the Schrodinger equation. One starts with the equa- 
tion 

fp 

(2.1) —^p{x, s) + q(x)^(x, s) — s'^^(x, s) seC 

dx'^ 

on the half line [0, -|-oo[. This equation has a unique solution, ip{x,s), satisfying 
the boundary condition 

(2.2) V'(0,s)=0, ^V(0,s) = l. 

For a potential q decreasing sufficiently at infinity there are also two solutions, the 
Jost solutions /(a;,±s), determined by the conditions 

(2.3) /(x,±s)eT^^^->™ 1. 
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For real fc, the Jost solutions /(x, k) and f{x, —k) correspond to incoming and 
outgoing waves and -0 is a linear combination 

(2.4) ^{x, k) = ^ {fix, k)M{k) - fix, -k)M{k)) 

such that M is the boundary value of the function M(s) = /(O, s), analytic in the 
upper half-plane. It follows that, as a; — > oo, 

Mx, k) = IMM sinlkx - r7(fc)) + o(l) 
k 

where ry(fc) = argM(fc) is the phase shift created by the potential q, with respect 
to the case g = 0. The problem of inverse scattering is, given the phase shift 
rj, to reconstruct the potential q. For phase shifts satisfying a growth condition at 
infinity, there exists a unique solution without bound states, which can be recovered 
by the methods developped in the 50's, see e.g. |12) and [9]. 

2.2. Discrete scattering. We now follow the exposition of Case and Geron- 
imo |8]. More information on orthogonal polynomials on the unit circle can be 
found e.g. in Szego |21j . Geronimus 1131, or the comprehensive treatise by Simon 

Let /i be a finite positive measure on the unit circle C/(l). Denote its moments 

e-^'^'dfiie) - Ck 

and (j)niz) the sequence of orthogonal polynomials obtained by the Gram-Schmidt 
procedure from 1, z, z'^ , . . ., normalized by 

(j)niz) = z" + ... 

This normalization will be more convenient for our purposes than the usual one, 
where the polynomials have i^-norm 1. Denote 



J \(j)nfdfj. and q;„ = 0„(O). 



The an are the Verblunsky coefficients of the system, cf |19j . They satisfy 

thus the sequence (T„ is determined by ctq = /i([/(l)) and the sequence a„. The 
orthogonal polynomials are subject to the recursion relation 

{(f>n+liz) = Z<l>n{z) + an+lZ'^(f>nil/z) 
_ 
0„+l(l/z) = Z "■(jjniz) + an+lZ Vn(l/z) 

In matrix form, putting 

(/>n(2) 



^iz,n) 
and 

Din) 

one can write 



0n(l/^) 



z an+iz' 

Un+lZ^" Z^^ 



^iz, n + 1) = Din)^iz, n) 
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Assuming that the potential a satisfies hm„ 
ost soluti 



,an — Q fast enough, one can 



construct the Jost solutions to equation (2.5) 

satisfying the boundary conditions 

lim|0+(z,n) - z"| = 

n 

lim|0_(z,n) - = 



n) 



(l)-{z,n) 



\z\>l 
\z\<l 
\z\>l 
\z\ < 1. 



lira \(f'+{z, n)\ = 

n 

lim 10- (z, n) \ — 

n 

The solutions are linearly independent, and 

(2.6) -^{z.n) ^ f+{z)-^+[z,n) + ,f_{z)^_{z,n) 

where 

/+(z) = /_(z) for |z| = l 

Equation (2.6 1 is the discrete analogue of (2.4 1. 

The function has analytic continuation inside the unit disc, it is related to 
the original measure on the unit circle by 

Me) = 

The phase r]{9) — arg(/_|_(e'^)) is the phase shift induced by the potential q;„. It 



can be recovered from the function 



on the unit circle. 



l/+(e'»)P 

One sees that reconstructing the phase shift from the potential a is equivalent 
to finding the orthogonality relation for the sec^uence of polynomials 0„, solution 
to the recursion equation (2.5 1. Conversely, the inverse scattering problem is the 
easy problem of finding the recursion coefficients for the orthogonal polynomials, 
knowing their orthogonality measure. In fact the connection between scattering 
theory and orthogonal polynomials has been known for a long time and is at the 
heart of the Gelfand-Levitan method (see I12| ). 

2.3. Comparing the discrete and continuous cases. Given a solution A 
of the Riccatti equation 

Q^A' + A^ 



the equation (2.1 ), for s = k real, is equivalent to the Dirac system 
(2.7) 



dy/dx 
-dz/dx 



= Ay + fcz, 
— Az + ky. 



The Dirac system (2.7) can then be put in the form 

AWe^'^'-", 



(2.8) 
with 

Similarly, putting 



dW/dx 
dW/dx = AWe-^^^"". 

W ^{y + iz)e'^^ 
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the system (2.5) becomes 
(2.9) 



a„+ie w„ 



It is now clear that (2.9 1 is a discretized version of (2.8). 

2.4. Recursion relations and the Caratheodory function. We recall a 
few properties of the orthogonal polynomials which will be needed in the sequel. 
The function 



F{z) 



-d^i{6) 



e'" - z 

is called the Caratheodory function of /i, it is holomorphic in C* \ and has 

the expansion 



z < 1 



Fiz) = co + 2Er=iCfe^ 
F{z) = -co-2S^^,c_fcZ-^- 1^1 >1 



Let 



and introduce the associated polynomials 

' e''^ + z 



■00 (^) 



= 1 



= 1 



+ z 



{z^Me'^)'K{z))dm 



and the functions 

e„(z) = i^n{z) + F{z)(l)n{z) 



e"" + z 
e*^ - z 



then, for n > 1, 

(2.10) e„(z) = 2(7„z" + C'(z"+^) z^O 

(2.11) en(2;) = 2(T„a„+i2:"i + 0(z"2) 
(2.12) 
(2.13) 



(t>{e^'>)z"'dij{0) 



oo 



eliz) = 2(T„a„+iz"+i + 0(z"+2) 
e*(z) = 2(T„ + 0(z^"'") z oo 



In other words, the fractions — ^ and ^ are reduced continued fractions for F, 
near oo and respectively. Furthermore, 

en+i{z) = ze„(z) - a„+ie* (z) 

= e,*(z) -a„+i2:e„(z). 
The recursion relations can be put in the form 

(t>7l{z) 



(2.14) 
(2.15) 



w(z) 



Yn{z) 

r„+i(z) = s„(z)r„(z) 
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with 
(2.16) 



za„+i 1 



There are also discrete Wronskian identities, derived from the recursion rela- 
tions, namely 

(2.17) 0„+i(z)V'„(z) - V'«+i(2)'/'n(2) = (/>„+i(z)e„(z) - e„+i(z)0„(z) 

2n 

(2.18) 0;+i(z)V';(2;) -<,+i(z)'/',*(2) </>;+i(z)e;(2;) - e;+i(z)0;(z) 

= 2a„+iCT„z"+i 

(2.19) 0„(z)V'*(2;) + V'n(^)'/'r*(^) = K{z)€a{z) ^ tn[z)(i>*^{z) 

= Ig^z^. 

3. Nonlinear relations for the Verblunsky coefficients 

3.1. Laguerre's method. We now adapt the method of Laguerre |16j and 
Freud |10j for finding nonlinear relations between the Verblunsky coefficients. The 
computations that follow are inspired by Forrester and Witte who treat the 
case of difference equations of the form wi^z + 1) = p{^z)w{^z). Another derivation 
for the (/-recurrence relations can be found in the paper |14| by Ismail and Witte, 
however we shall need them in the form below, which is slightly diflterent from the 
results of [14j . so we give a detailed derivation here. Let us assume that the measure 
\x has the form 

d\x{d) = w{e^)dQ 
with w holomorphic in C*, satisfying 
(3.1) w{qz) = p{z)w{z) 

where q g]0, 1[, and p is a rational function, 

for some polynomials V, W . The case of interest is when 



(3.3) ^(6*") 
q e]0, 1[, a,& e C*, thus 

satisfies 



(ae^','Z)c 



(6e^«,<?)c 



{ae^^ ,q)oo{ae '^,g)c 
{be-^,q)oo{he-'',q)o 



{az,q)oo{az \(j)o 

w[z) — 



{hz,q)oo{hz 1,5) 

C 



(qz — o){az — 1) 

and we may take 

(3.4) V{z) = {qz-a){bz-l) 

(3.5) W{z) = {qz-b){l-az). 
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For the moment we shall make the computations with general polynomials V 
and W, and come back to the specific values (3.4), (3.51 in the last section. 
Equation (3.1) implies that the Caratheodory function satisfies 



(3.6) 



W{z)F{qz) = V{z)F{z) + U{z) 



with U a polynomial of degree at most D — max(deg VF, deg y). Indeed, take the 
partial fraction expansion 



piz)^6iz) + J2-^ 



3 

z — a," 



S polynomial of degree less than (deg^ — deg W^). 



with W{z) = Cl\j{z ~aj), then 



by a change of integration contour 



p{z)F{z)+J^^JpiC)-p{z))wiC)§ 



^ piz)F{z) + J 5(C,z) + E 



j z~aj C- 



_C+£ 



~S(C,z)^{C + zf-i^^ 



dc 



{Viz)F{z) + U{z))/Wiz) 



where [/ is a polynomial of degree at most D. We shall now use (3.6) and the 
characterization of the quotients ipn/4'n and V'n/^^n as reduced continued fraction 
expansions. Let us start from 

e„(z) = Ipniz) + F{z)(l>n{z) 

and 



eniqz) - Mqz) + i^Fiz)+^^'^ 



(qz) by (3.61 



'Wiz) ' ' Wiz) 
and eliminate F{z) between these two equations, to get 
W{z)en{qz)(l)n{z) - V{z)eniz)(t>n{qz) = 

W{z)lpn{qz)(l)n{z) - V{z)^pniz)(t>n{qz) + ?7(z)(/i„ (z) (/)„ (gz) . 



The right hand side is a polynomial in z and, according to ( 2.10 )-( 2.13), the left 
hand side is Olz^) near and 0(z^+"-i) near oo, therefore there exists a polyno- 
mial Qn{z) of degree less than D — 1, such that 



(3.7) 



2a„+icr„z"e„(z) = W {z)en{qz)(t>n{z) - V{z)en{z)(t>n{qz). 



In case a„+i = 0, one has 0„+i(z) = z(/)„(z) and 4>n^i{z) — (j^niz), and the sought 
relation for (f)n+i{qz) can be obtained from that of (j)n{qz). From 



2a„+icr„2;"e„(z) = W{z)'ipri{qz)(j)n{z) - V{z)ipn{z)(j)n{qz) + U{z)(pn{z)(l>n{qz) 
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we get, according to (2.17), 

[(f>n+l{z)'ljjn{z) - ^pn+l{z)(t)„{z)]Qn{z) 



Wiz)ijniqz)(j>n{z) - V{z)'il;n{z)(l)n{qz) 

+ U{z)(f)niz)(f)n{qz) 



V'„(2:)((/)„+i(z)e„(z)+V(z)0„(9z)) = (l)n{z){W{z)i^n{qz)+U{z)(j)Jqz)+il;n+iQn{z)). 

This expression is a polynomial, common multiple of and ipn which have no 
common zero (see (2.17)), so it can be put in the form 

^n(z)(l)n{z)llJn{z) 

where r2„ is a polynomial of degree at most D. Thus Qn satisfies 

(3.8) r2„(z)0„(z) = (t)n+l{z)en{z) +V(z)(f>n{qz) 

(3.9) n,Az)ijn{z) = W{z)MQ^) + U{z)cf>n{qz)+^n+l{z)en{z) 



One obtains ri„ by multiplying (3.9) by 4>n+i, substracting (3.8) multiplied by 
ipn+i, and using (2.17), thus 

(3.10) 2a„+i(T„z"ri„(z) = W{z)en{qz)(t>n+iiz) - y(z)e„+i(z)(/)„(gz) 

jives 

V{z)(f)niqz) = ri„(z)0„(z) - (/)„+i(z)e„(z) 



Equation (3.8) gives 
(3.11) 



or, using (2.5), 
(3.12) 



V{z)4>n{qz) = i^niz) - ze„(z))0„(z) - a„+i e„ (z)(/>* (z) 



On multiplying (3.81 by F{z) and adding (3.9) 



Q.niz)e„{z) = V{z)en{qz) 



w{z) 
w{qz) 



+ e„+i(z)8„(z) 



therefore 

Finally 
(3.13) 
and 
(3.14) 



= vizf-^ + ^^e„(z). 



w{z) 



w{qz) w{z) 



v{zf-^ ^ n^izf44 - '^e„(z) 



w{qz) 



w{z) w{z) 



V{z) 



- (S2„(z) - ze„(z))— — + a„+iB„(z)- 



w(qz) W[Z) W[Z) 

Equations for 0* and e* are obtained by similar arguments. One finds polyno- 
mials 6*1 (z) and f^^(z) with the same bounds on their degrees as 6„(z) and f2n(z), 
such that 



(3.15) v{z)(j)*„{qz) = (f7;(z) - e;(z))(/);(z) - za„+ie;(z)(/)„(z) 

and 



(3.16) 



vizfS^ = mz) ~ e:(z))^ + za„+ie:(z)-^"(") 



w{qz) 



w{z) 



w{z) 
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Proposition 3.1. Let 



(3.17) 



then, withYn defined by {2.14.), the relations [3.12), {3.14), (3.15), (3.16) can 
he put in matrix form as 



(3.18) 



Yniqz) 



Viz 



-Ar,{z)Yr,{z) 



Compatibility relations between (2.15) and (3.18) imply that 



(3.19) 



An+i[z)Bn{z) = Bn{qz)An{z). 



The coefficients of the polynomials 8„, 0* , ri„, $7* in (3.171 can be computed 
from the coefficients of the polynomials 0„. For example, with V and W given by 



(3.4) and (3.5) one finds, putting 



9„(z) = Az + M 



(3.20) 2a„+iz"e„(z) = W{z)e^{qz)4>Jz) - V{z)e^{z)<t>^{qz) 

and comparing the two expressions near zero, 

2a„+iz> = 2a„z"(6g"-a) 
Comparing at infinity gives A, and finally 



(3.21) 



e„(z) = (a - hq''+')z + {-a + fog")- 



an+i 



Similar reasoning yield the other polynomials in terms of the Verblunsky coeffi- 
cients. Using 

0„(z) = z" + /3„z"^i + ...+Q„ 



with 



3„ = 

j=i 



One finds 
(3.22) 



r2„(z) - ze„(z) = + tz + bq'' 

where t is a complicated expression in terms of a„, q;„+i and /3„. Similarly 



(3.23) 
(3.24) 



Q*Jz) = [aq - 6g"+l)^^z + - d 

ri* (z) — 6Jj(z) — aqz"^ + t*z + d 



Equation (3.19) implies nonlinear recursion relations among these coefficients. 
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3.2. The isomonodromy deformation. Rather than do the tedious compu- 



tations here, we shall, in the next section, start from the relation (3.191 and derive 



everything from its qualitative features. This will allow us to make the connection 
with Sakai's theory of discrete Painleve equations. Before this we first make a few 



remarks. Equation (3.191 has the form 
(3.25) 



where 



and 



A{z) = B{qz)A{z)B{zy 



A{z) = Adz'^ + ... + Aq, A{z) = Adz"^ + . . . + io 



B{z) = B 



z 
1 



for some constant matrix B. Equation (3.25) is an isomonodromy transformation 



of a q-difference equation, as studied by Birkhoff [2], and recently by Borodin 
[3]. It is closely related to the equation appearing in |15| . Such isomonodromy 
transformations, for difference equations instead of g-difference equations, have 
been studied in depth by Arinkin and Borodin in [1 , using geometric methods. 
These methods, involving the definition of "g-connections" should also apply to our 
situation, although we do not use them here. 

In our case, the matrices Ad, Ad are lower triangular matrices, and ^0,^0 a-re 
upper triangular. Let ki,K2 be the eigenvalues of Ad, and 0i,92 those of Aq, so 
that 



«2, 



Ad 



Kl 



An = 







Consider the matrix 



then 



Ad = 



Aiz) 



qui 




qz 




H2 



Aiz) 



Ao^ 



qOi 



= A2Z^ + Aiz + Ao 



detA{z) = qdctA{z) 



We can now conjugate A{z) by a constant matrix B in order to get the matrix 
A{z) with 

'qui \ 4 ^ * 

* K2 r ° I 02, 



Ad 



The matrix B is determined up to left multiplication by a diagonal invertible 
matrix. One can choose its coefficients to be expressed as rational functions of the 
coefficients of A{z). Furthermore equation ( 3.25| can be used in the same way to 
deduce A{z) from A{z). It follows that equation (3.25) defines a birational transfor- 



mation between two spaces of matrices with polynomial coefficients. We shall give 
an explicit formula for this transformation in the next section. The recursion equa- 
tions for Verblunsky coefficients are given by iterating this transformation. In the 
next section we shall investigate in details the case where d = 2, which corresponds 



to the case of the Verblunsky coefficients for the g-Gamma weight (3.3) 
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4. A birational transformation 

4.1. The space -'i^Ki,K2,ei,6i2,ci,c2,c3,c4- We now embark on the analysis of equa- 



tion (3.25 1, in the case of degree 2 polynomial matrices. Consider a two-dimensional 
vector space iJ, and let 

A{z) ^Aq + Aiz + A2z^ 

be an endomorphism of H which depends polynomially on the variable z, with 
^0,^1,^2, elements of End(i7). 

Definition 4.1. We let 

■^Ki,K2,0i,02,Ci,C2 .03,04, 

denote the space of all such A{z) satisfying 

{ki,K2} ^ spA2, {61,62} ^spAo 

and 

(4.1) detA{z) — KiK2{z — ci){z — C2){z — C3)(z — C4) 

for some fixed complex constants ki, K2, 61, 62,01, ... , C4, satisfying the equation 

(4.2) K1K2C1C2C3C4 = 6162 
(which comes from considering detA{0)). We let 

^Kl,K2.9i,6'2,Ci,C2,C3,C4 

denote its quotient by the action of GL{H) by conjugation. These spaces are com- 
plex manifolds. 

In the sequel ci , . . . , C4 will be fixed so we will omit them from the notation 
most of the time. 

We consider the case of generic parameters, in particular ni ^ H2,6i ^ 62, and 
^0,^2 are semisimple. We denote Ei,E2 the eigenspaces of A2, and Fi,F2 those 
of The set of A{z) such that E2 7^ Fi is a dense open set il, invariant by 
conjugation, and for each element of this open set we can chose a basis (u, v) with 
u £ Fi,v Cz E2 so that, in this basis, Aq and A2 have matrices of the form 

\* K2J \0 62 

As in the preceding section, let g be a nonzero complex number and consider the 
matrix 

M^) = J) M^) ('q J) = A2z' + Alz + Ao 



then 



A2=(7 *], Ao^r' 1), detAiz)^qdetA{z) 



72 



therefore the matrix A{z) belongs to the space Xqi^^^i^^^gg^^g^. Thus A{z) A{z) de- 
fines a map from a dense open subset of Aki,k2, 61,62 to A^Ki,K2,gei,e2 ■ Furthermore, 
as explained also at the end of the preceding section, this map can be inverted, thus 
defining a bijection between dense open subsets of Xf^^^^f^^, 01,62 and A'gKi,K2,9ei,e2- 
Clearly the map passes to the quotients and induces a bijection between dense open 
subsets of A„j ,.;2, 61,62 ^^'^ ^9^1,^2, 961,62; which is a birational map. 
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Definition 4.2. We call 4> : X^^^^^^^fi^ -> ^9Ki,K2,9ei,e2 
structed above. 



the birational map con- 



In the next section, we shall introduce a parametrization of the spaces in order 
to give an explicit formula for this birational transformation. 

4.2. Jimbo and Sakai parametrization. We will use the following parametriza- 
tion of the space ^Ki.K2,ei,02 1 Stitev Jimbo and Sakai [15j 



(4.3) 



A{z) = 



Ki{{z - y){z - a) + zi) z-y 

K2{{Z - y){z -S) + Z2) 



■4iz + /3) 

Thus y is the unique root of the polynomial which is the 12 coefficient of A, and 
KiZi, K2Z2 are the values of the 11, resp. 22 coefficient of A{z) aX z = y. One has 

ya + zi = Oi/ki, ya + Z2 = ^2 /k2 

The quantities y, zi, Z2 satisfy 

(4.4) Z1Z2 = (2/ - ci){y - C2){y - c^){y - C4) 

where the equality follows from considering the determinant of A[y). Let us put 
S,^{y- ci)(y - C2)/kiZi = K2Z2/{y - C3){y - C4). 



The quantities /3 and 7 are determined as rational functions of y, ^ using (4.1 ) 



It follows that the variables (y, ^) in C* x C* parametrize a dense open subset 



Ki,K,2,Oi.,H2 ■ 



Using this parametrization, the space X^^^ .^^ 



is identified with 



of X, 

X P-'^ blown up at the seven points 

(Ci,0), (C2,0), (C3,00), (C4,0o) 
(O,CiC2/0l), (O,CiC2/02), (00,1/ki) 

minus the strict transforms of the four lines y — 0,y — (X),£, — 0,^ — 00, see Fig. 2 
below. 



f = CXD 



(0, cva/e. 



(O,ClC2/02) 




(00, 1/k.i) 



(Cl,0) (C2,0) 

y = y = 00 

Fig. 2. The space 



The curves corresponding to F2 = Ei and F2 = E2 are respectively the strict 
transforms of the points (0,CiC2/6'2) and (0, CiC2/6'i). 
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4.3. Explicit birational expression for (p. We now describe the map (j) ^ 



a birational transformation, using tlie Jimbo-Sakai coordinates (4.3 1. Let A{z) be 
as in ( |4.3[ ) and consider the matrix 

^ ( qni{{z - y){z - a) + zi) qz{z - y) 

\ (72 + /9) K2{{z ~ y){z - 5) + Z2) 

Then A{z) belongs to the space XqKi,K2,g6»i,6i2 • In order to find its Jimbo-Sakai 
coordinates, we have to use the basis formed by the spaces E2 and Fi. It follows 
that 

A{z) = B-^A{z)B 

with 



B = 



H ~ O2 ~q 

/3 QKi - K2 



on the condition that the matrix B so defined is invertible. In order to go from A 
to A one proceeds in the opposite way. Start from 



A^ 



qKi{{z ~ y){z - d) + zi) z-y 

z{jz + P) K2{{z - y){z - S) + Z2) 



conjugate it by a constant matrix, so that its 21 coefficient has degree 1 and its 12 
coefhcient is zx(a polynomial with degree 1). For this use eigenvectors of A2 and 
Ao, thus A{z) = C-^A{z)C with 



C 



qni - K2 -y 
97 qOi - 6*2 



We will not compute the inverse transformation. 

Let us now compute an explicit formula for (j), using the coordinates {y,0- We 
start from the matrix 

/kiz'^ + az + 01 z-y \ 

^^""'-y z(7Z + /3) K2Z^+bz + 02j 



with 



(4.5) a = -KiyH 

y 

(4.6) b = -K2yH 

y y 

and 



(4.7) 0ib + 02a + (3y = -kiK20-3 

(4.8) K16 + K2a — 7 ~ — KlK20'l 

where the Ci are the elementary symmetric functions in the c^. Then 

-T/N „/ 1 f qK-i z^ + dz + q0-\ w(z — y) 

4.9) A(z) = Biqz)A(z)B(z)-^ = ^ , 3^ 2 I 

^ ^ ^ ^ ; V ; V ; ^ w^^z{jz + j3) K2Z^ + bz + 
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B{z) = 

Biz)-' = 
A ^ 



with 

^z{qKi - K2) 

-z(3 

A computation using (4.9) gives wy = q, and 



{qni 



^R/A, 



y 



(71 - 

O2)- 

qA 
R 



qui 



Q 

K2 



with 



iy^A = CyHqOi-e2){qK,-K2)+qiy^P 
= (,y^{q0i - 02){qKi - K2) 



and 



+9y^(-KiK2cr3 - 6*16 - 6*20) usingpT? I 
-qKiK2Bi(^{y - C3){y - C4) 
+^{{q^Kidi + K202)y'^ - qKiK2(T3y + 2q0i02) 
-q02{y - ci)(y - C2) using K^, (|46]), 



^yiiqi^i - H2){~aq^ + bq - q{qKi - K2)y) - q^l) 
iy(-q^Kia ~ qn2b - q{qni - Ka^y - q^niK2(7i) by (|48| 
q^niim^y^ - (y - ci)(y - cs) + 

+qi^2{i^2iy^ - KiK2£.^{y - cz){y - C4) + 6*20 

-q{qni - K2fy'^i - q'^KiK2<Ji£,y by ( |X5| ), ( |4!6l ) 

-qninl^'^iy ~ C3){y - C4) 

+^(2g^KiK22/^ - q^KiK2criy + q{q^KiOi + K2O2)) 
-q^i^iiy - ci)(y - C2). 



Therefore 
(4.10) 



y 

with 
S 



yT 



qKlK2<J3y + 2q0i02 



-qKiK20i(.'^{y - C3){y - C4) 

Hiiq^i^idi + K202)y^ 

-q02{y - ci){y - C2) 
T = -K,inl^^{y ~ C3){y - Ci) 

+^{2qKiK2y^ - qniK2(Tiy + {q^Ki0i + K202)) 
-q^Kiiy - ci){y - C2). 
In order to find ^ we proceed as follows. From the definition one has 

{y - ci)(y - C2) 



qKip +dy + q0i 



where a can be computed from (4.9). Plugging in the values of a,b,(3,j, and y 
given by (4.5), (4.6), (4.7) (4.8), (4.10) gives an unwieldy expression as a rational 
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fraction in the variables y,^. One can simplify this expression using the following 



observation. Computing the determinant detA{y) = qdetA{y) from (4.9) gives the 
identity 

(4.11) qKiK2{y - ci){y - C2){y - C3){y - Ci) {qKiy"^ + ay + qei){K2y'^ + by + 02) 

where y, a, b are rational expressions in y, ^. Of course this identity is compatible 
with the alternative expression 

~ _ K2y^ + by + e2 



qKiK2{y - C3){y - a)' 



Replacing y, a, b in (4.11) by their values and reducing the denominators, we get 
a polynomial identity between the numerators. This leads us to expect that the 
numerator of each expression y — Ci factorizes into a product ot two terms, which 
divide respectively the numerators of quiy^ + dy + q9i and K2y^ + 6y + 02- Since the 
degree of the numerator of y — Cj in the variable ^ is two, we expect a factorization 
into two terms of degree 1 in i^. Indeed a direct computation confirms that 

S-ayT = KiK2^'^{ciyK2 - qOi){y - C3){y ~ C4) 

+C [-2qciKiK2y^ + {q^KiOi + K262 + qciKiK2(Ji)y'^] 
+C [{-qKiK2(T3 + ci{q^ni0i + K292))y + 2(76'i6'2] 
+q{qciKiy - 6'2)(y - ci)(y - C2) 
= {^{ciyK2 - qOi) - qci{y - C2)) ^ 

x(KiK2C(y - C3){y - Ci) {qcimy - 92){y - ci)) 

Cl 

and similar expressions for the other Ci, e.g. 

S-C3yT = KiK2^'^{c3yK2 ~ q9i){y ~ C3){y - C4) 

+^ [-2qc3KiK2y^ + {q^ni9i + K292 + qc3KiH2<yi)y'^] 
+C [{-qi^iii2<y3. + C3{q^Ki9i + K2^2))y + 2g6'i6'2] 

+q{qc3Kiy - 92){y - ci){y - 02) 

= {£.{y-c4) {qcsKiy - 62)) -x 

KIK2C3 

x(KiK2C(y - C3){c3yK2 - q9i) - KiK2C3{y - ci){y - C2)) 

One checks then that one of these factors divides the numerator of quiy^ +ay + q9i. 
This yields the simpler expression 

.4^2. C1C2 iay~£t.)~^^y~c2))iay-£k)-^.(y-'^'^)) 

n,9,^ay - c,) - i^{y - ^j){ay - C3) - i^{y - j:^jy 



Proposition 4.3. Formulas {JjJ^ o.nd {4- 12) express (j) as a birational transfor- 
mation (j) : {y,C) 

4.4. Blowing up the spaces. The map 0, whose expression as a birational 
transformation we obtained in the preceding section, can be extended to a regular 
map between two surfaces, by blowing up the spaces ^Ki,K2,6ii,e2 ^^'i Xqiii,K2,q9ifi-2- 
For this we first extend the map to the set of endomorphisms such that E2 = 
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Fi. As noted before, this set corresponds to the blown up point with coordinates 
(0, Ci 02/^2)- Take a matrix in a neighbourhood of E2 = Fi in the form 



Kiz'^ + az + 62 z — y 

z{jz + i3) K2z'^ + bz + 



with y small and ^ — ciC2Ki/02 = = [y — ci){y ^ C2)/zi, or 

Kiy - Ki(ci + C2) - \62 
a = -Kiy H r — ■ 

Ay + CiC2Ki/e2 

We can conjugate the matrix to put it in the form 

Kiz'^ + az + 01 z — y 



with 



then putting 



Z{^Z + /3) K2Z^ + hz + i 
_ 02 — 01 

a = a H 

y 

b = 1+^-^ 
y 

^ - , (ki - K2){01 - 02) 

7 = cri + KiO + K2aH 

y 

^ a3 + 0ib + 02a {01-02)^ 

y y^ 

A = {q0i - 02){qKi - K2) + ql3 
A 



y = 



{qki - K2)i-qa - y{qKi - K2)) + (-97 + (^^i - «2)&) 
we obtain 

- _ (g/ti - K2){q0i - 02)qa - gyPiqui - K2) + g7(g^i - ^2) + qi3b 

A 

We keep only the main terms as ?/ — )• to get 

A ~ -gi^ 

" ^ (g^Ki-re2)!/ 



a ~ 



therefore 



-2 I ~- - VCl ~ ^2 ~2 I 

KiZi ^ qKiy +ya^ ^22/7^ ^ ^ H2y q 

{q^Ki - K,2)y 



or 

f = 9/K2 

Thus we see that the map can be extended if wc blow up the final space at (00, q/ K2). 
A similar computation shows that the inverse map also extends to the curve E2 = Fi 
if we blow up the initial space at the point (00, q/K2)- 

Proposition 4.4. Let 5^K;i,K2,ei,02 denote x P^ blown up at the eight points 

(ci,0), (C2,0), (c3,oo), (c4,oo) 

(0,CiC2/6'i), iO,CiC2/02), (00,1/Ki), (00,g/K2) 

minus the strict transforms of the four lines y = 0,oo, ^ = 0,oo. then the map 
extends to a regular isomorphism between the spaces Y^^^^^^g^^g^ and l^rei,K2,9ei,02- 
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^ = 00 



(0,ClC2/9l) 




(00, 1/ftl) 



{co,q/K2) 



{Cl,0) (C2,0) 

y = y = oo 

Fig. 3. The space ¥^^^^^0^92 
5. Discrete Painleve equations 

5.1. Sakai's theory. We recall here the basic facts about the relation between 
surface theory and discrete Painleve equations, according to Sakai |18j . One starts 
with a complex surface X, obtained by blowing up P^(C) at nine points xi, . . . ,xg, 
thus X comes with a morphism tt : X — > P^(C). Degenerate configurations, where 
some of the points Xi are infinitely close, are allowed. The Picard group of such a 
surface is a free abelian group of rank 10, generated by £q = Tr~^{L), the inverse 
image of a generic hnc in P^(C), and Si = ■K~^{xi), i = 1, . . . ,9. The intersection 
form on Pic{X) satisfies 

'^^o-'^^o = 1, = -1, j = 1, . . . , 9, £i.£j =0, i ^ j 

The canonical divisor class is 

K:x = ~6^ ~3£o +£i + ... + £9 

and its orthogonal S-^ in Pic(X) is a rank 9 lattice, generated by 

which form the simple roots in an affine root system of type Eg^"^ 
Dynkin diagram 

£0 ~ £1 — £2 — £3 
Q 



according to the 



£1 — £2 £2 ~ £3 £3 — £i £i — £5 £5 — ^6 ^6 ~ £7 £7 ~ ^8 

Fig. 1 Dynkin diagram of the affine root system E^^^^ inside Pic(X) 



— o 

^^8 — £9 



The surface X is called a generalized Halphen surface by Sakai, if the linear 
system | —ICx \ has a unique divisor, whose irreducible components are orthogonal to 
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ICx- Let i G / be these components, then they generate an afBnc subroot system 
i? in (5-'- C Pic(X), such that its orthogonal is also an affine subroot system C 
(J-*-. These surfaces are classified according to their corresponding root systems. Let 
X and X be two such surfaces, corresponding to the same root system R, where the 
irreducible components of the canonical divisor are labelled « S /, resp. Vi^i € I, 
then Sakai's main result is that (up to some little twists, see |18j . Theorem 25 for a 
precise statement), every isometry of the Picard groups 4> : Pic(X) — > Pic{X) which 
preserves the components of the anticanonical divisor, (t){{'Di,i e /)) = G /) 

comes from a unique isomorphism of surfaces (jj : X X. Such isometrics are 
given by elements of the extended Weyl group of the dual root system (i.e. the 
group generated by the Weyl group of R-^ and the automorphisms of the Dynkin 
diagram). When viewed at the level of P^(C), the isomorphisms (f> : X ^ X give 
rise to birational transformations P^(C) P^(C). The discrete Painleve equations 
correspond to the translations in the affine Weyl group W{R-^). On the other hand, 
transformations in the extended Weyl group of R give Backlund transformations. 
As an example, Sakai treats the case of the birational tranformation 

P ■■ if, a) ^ (7,5) 

determined by 

_ ^ hb^ if-b,){f-be) 
^ g if-brKf^bs) 

- ^ ^7^8 {g - qbi){g - qb2) 

f {g-h){g-bi) 

where 6i , . . . , 5g are some complex parameters and 

_ b-jbjb^bQ 

He shows that this transformation corresponds to an isomorphism between A^^^ 
surfaces, which arises from a translation in an affine Weyl group d'^^\ Here Ag^-* 
and Dg"'^'' are two orthogonal affine root systems in e''^\ Remark that this trans- 
formation is actually considered as a birational tranformation of P^(C) x P^(C). 
Since P^(C) blown up at 9 points is isomorphic to P^(C) x P^(C) blown up at 
8 points the two pictures are equivalent, but it is sometimes easier to do formal 
reasoning on P^(C) and to perform actual computations in P^(C) x P^(C). 

As explained in |15j . the transformation p can be considered as a step in a 
nonlinear difference equation 

{fn+l, gn+l) — P{fn,gn) 

and this difference equation leads, by taking a suitable continuous limit, to a rational 
differential system 

df/dx p{x,f,g) 
dg/dx = q{x,f,g) 

which is equivalent to a Hamiltonian version of the Painleve VI equation. The 
degenerescence scheme of discrete Painleve equations corresponds to the inclusion 
of the root systems in E^^\ The degenerescence to Painleve differential equations 
is also discussed by Sakai. 
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5.2. According to the classification of Sakai, the spaces yKi,K2,Si,S2 *3x- 
amples of Ag^^ surfaces. We wih now identify the map (p between yKi.K2, 61,62 
and yqKi,K2,g6»i,e2 have constructed as one of the discrete Painleve maps. For 
this we will use intersection theory for surfaces obtained as blow-ups. Denote by 
J^i,i — 1, . . . ,8 the divisors which are inverse images, respectively, of the points 

(Cl,0), (C2,0), (c3,oo), (c4,oo), 
(0,CiC2/6'i), {O,CiC2/02), (00,1/ki), (00,q/K2) 

in the blow ups. Similarly we denote hy J- i,i = 1,...,8 the ones in the space 

^qiil ,K2 ,q6i ,82 ■ 

It will be convenient to identify, as in Sakai [181, the surfaces with blown 
up at 9 points, called Xi. This identification is as follows. Consider x P^ with 
coordinates denoted by {f,g), blown up at the eight points 

(0,61), (0,62), (00,63), (00,64), 

(^5,0), (65,0), (67,00), (68,00) 

We make the explicit identification 

Cl = 61 C1C2/61 = 65 

/r- IN C2 = 62 C1C2/O2 =65 , ^ 

^^■'^ C3 = 63 1M=67 ^ = ^ 

a = 64 q/ K2 = 68 

between our space and Sakai's. Consider also P^ blown up at the eight points 

xi = (fli : : 1), a::2 = (02 : : 1), 2:3 = (03 : : 1), 0:4 = (0 : 1 : 04), 
0:5 = (0 : 1 : 05), Xq = {1 : ae : 0), X7 = (1 : 07 : 0), X8 = (0 : 1 : 0) 

further blown up on 7r^^(0 : 1 : 0) at 

2^9 = ( -) = (0,08) 
\z X / 

The two spaces are identified according to the transformation 
(5.2) 

61 62 63 64 „ \ /l/a2 1/03 1/ai as V z_ 

65 66 67 68' ''^y \ ae 07 -l/a4ai -l/asoi ' a: - aiz ' x 

See pages 173 and 215 in |18j . The correspondance between divisors is 

£t = J-S, £g — J-4, £0 — £1 — £& = J'3 

£1 is the line y = a^, and £§ ~ £q the line ^ = 00 . 
The components of the anticanonical cycle are 

2^0 = ^8 ^ £9 T^i = £0 — £& — £7 — £& 

T^2 — £0 ^ £1 ^ £2 ^ £3 T^3 — £0 ^ £4 ^ £5 ~ £& 
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which form a root system of type A'^\ with Dynkin diagram 




The orthogonal is a root system of type dI^^\ with generators 



(5.4) 



ao 
as 



So 



£i-£s 



£91 
£&i 



forming the Dynkin diagram 



ai = £2 -£3, a2 = £i-£2 

0^4 = £q—£7, 0^5 = f 4 — f 5 




We denote by £ and F the objects corresponding to the target space of the map (j). 
We saw that maps J^g to J-g, and J-^ to 7^5. We will compute the matrix of the 
map induced by between the Picard groups. 

5.3. Let us start by computing the image of the curve Ti = £2- This curve 
consists of matrices of the form 



A{z) 



Kiz^ + az + 61 w{z — ci) 

'W^^z{"fZ + (3) K2{z - Cl){z ~ 92/C1K2) 



where a is a local parameter along the curve and w an irrelevant constant which 
can be removed by conjugation with a diagonal matrix but is convenient for com- 
putations (this is the "gauge freedom" in [15 J. The parameters 7 and /3 are afhne 
functions of a. Recall that 



A = {qui - K2){qOi - 02 



qf3 



then 



y = 



qa{qKi - K2) - {qui - K2)(-ci - Oi/aiKi) + ci(gKi - K2Y + ^{qui - K2) 



The critical value to be investigated are y = 0, Ci, 00. The case y ~ corresponds 
to El = F2 so there is an intersection with J-q — £5. In the case y — ci, the matrix 

A{ci) has the form | £ ? |, therefore 



T 



i(ci) 



f{qei-e2)w -qW 

{{q0i-02)Z -qZ 
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{qr + qKi - K2)((qKi - k.2)G - 
[qr + qni - K2){I3G + qOi - 



qH) 
'2) 



with W — {qui — K2)qS + I3T, thus y — ci ii W — 0, and then zi = 0, therefore the 
hxiage does not mterscct J-i = £2- In the case y — Ci, with i = 1,2, 3, the matrix 
A{ci) has rank 1 so 

, frG 

for some constants H and G (independent of a), and 
AMi - e2)r - /3)((gKi - n2)G - qH) 

\ iq9, - e^y - i3)ii3G + qe^ - 62) 

The 12 coefficient vanishes only if qr + qKi — K2 = 0, and then Z2 — 0, therefore 
the image intersects J-i ii i = 2 but not i = 3, 4, and intersects £1. Finally for some 
value of a the equation y — 00 holds and the image intersects J^y. We conclude 
from these considerations that, in the Picard groups of the surfaces, 

(t){£2) = 2£q - £1 - £3 - £5 ~ £q ~ £s. 

Similar tedious computations, which we omit, lead finally to the matrix induced 
by the map (jj between the Picard groups, which is, from the basis {£i)i=o,...,9 to 
i£i)i=o,...,9 given by 



/ 6 


2 


2 


2 


3 








3 


2 


1\ 


-2 





-1 


-1 


-1 








-1 


-1 





-2 


-1 





-1 


-1 








-1 


-1 





-2 


-1 


-1 





-1 








-1 


-1 























1 











-3 


-1 


-1 


-1 


-2 








-1 


-1 


-1 


-3 


-1 


-1 


-1 


-1 








-2 


-1 


-1 

















1 














-2 


-1 


-1 


-1 


-1 








-1 



















-1 








-1 








One checks that this transformation preserves the components Vi (cf (5.3|) 



0(I?o) = 



V2 

T>3 



</.(I?3) = A 



while its action on the orthogonal root system ( |5.4[ ) is a translation in the affine 
Weyl group, namely 



0(qo) = ao, 0(ai)=ai, 0(^2) = "2, 
(j){a3) = as, (/)(a4) = a4 - (5, ^(as) = 0:5 + 5 

with 

S = ao + ai + 2a2 + 2as + + 

This translation has the following reduced decomposition, where the Wi are the 
simple reflexions associated with the roots a^, and a is the automorphism of the 
Dynkin diagram such that cr(ao) = ai, o'(q!4) = 0:5. 

(5.5) if) = aW4W3W2WoWiW2W3W4 

We summarize the results of this section in the following. 
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Proposition 5.1. The transformation cj) induces, on the Picard group of the sur- 
faces, a translation whose reduced decomposition is given by 

if) — aWiW^W2WQWiW2W^W4 

Remark 5.2. The discrete Painleve equation of |15) corresponds to the translation 
ip given by 

'4){ai) = Ui, i = 0, 1,4, 5 
■0(a2) = a2 - (5, ^{a3) = a3 + S 
The translations and tp are not conjugate in the Weyl group. 



5.4. As a check, we will obtain formulas (4.101, (4.12) using (5.5) and the ex- 



plicit expressions for the birational transformations corresponding to the reflexions. 
After Sakai, we have the following formulas for the transformations corresponding 
to the roots 

/ as ai 02 as 

Wo : 

\ 04 05 flg 07 
'l/ai 1/as 02 03 
^0104 a^ai a^as a^as 

wi,W4,wc, are the transpositions of parameters 

: (02,03); W4: (04,05), W5: (05,07) 



■,x:y: z 

x{x — Oiz) : y{x — z/as) : z{x — aiz] 



Passing to x P^ with (5.2), W2 and W3 are 

%i 62 63 64 
.65 



W2 



b2 
be 



63 

bi 



W3 



bi 62 

'bi 62 ^63 
ybr be 65 
The automorphism a is the inversion 
^bi 62 63 64 



63 

67 ^8 

bi 64 

67 bs 



63 
67 



;/,.9 

J g-bi 

;/,.9 



/-&7 



62 63 

be bs 



'Mbs/f, b^bi/g 



We can now compute </> induced by the translation (5.5). One gets 
6162^8 (/(.9 



/ 



65 



&8(.9-6i))(/(5 



65 ' 



b?.{g 



b^f {fig ~ 63) - 68(.9 " '^)(/(5 - ^4) - bsig t^) 



bs 



b,bj f-k 



9 f-bs 



-, f-l'8 



-b. 



This gives (4.10), (4.12), using the correspondance (5.1) and after simplifying the 



expression for g, noting that / — 65 divides the denominator, while f — bs divide 
the numerator. 
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5.5. Finally the results of the preceding section allow one to derive relations 
for the Verblunsky coefhcients. Indeed one has 

det An{z) = q'^V{z)W{z) 



thus one can take ci — a/q, ci — \/b C3 — b/q, C4 = 1/6. From (3.21) o ne see s 



that the coordinate y corresponds to ^ ""(a-6q"+i) • According to equations (3.221, 



(3.24), one finds that ni = bq"'^^, K2 = aq, 61 = 62 — a- We conclude 



that these expressions satisfy a discrete Painleve equation. 



6. Recovering the differential equation 

We shall now explain how to obtain a differential system as a limit of the 



discrete dynamical system given by (4.10), (4.12). For this we introduce a complex 



variable t, and consider the parameters [tni, K2,t6i, 62, €1,02,03,04). Consider the 
map (j) of Proposition (4.3) as a birational map {y{tq),£,{tq)) — (t){y{t), ^{t)). Using 
the parametrisation 

q=l-e, 
qKi = t{l + eKi), K2 = l + eK2 
g^i =<(l + eei), 92 = 1 + 692, 
Ci = 1 + eCi,i = 1, ... ,4 
y = 1 + eu, ^ = w 

Letting e go to zero gives a differential system in the variable t. We omit the details 
of the lengthy computation, and give the final result: 



t{l-t) 
t{l-t) 



du 
~di 

dv 



= tv{u - C3)(m - C4) - t{t - 1)u-1(m - Ci){u - C2) 

= v[2u + 2t{K2-Qi) + Ci+C2-K^-Qi] + 
2u{t - 1) + C2 + Ci - t{C3 + C4) + 
y-^ [2ut - t{l + K2) + 2(1 + 62) - C3 - C4] 



Finally, we come back to the original problem about the Fourier transform and 
the scattering problem for phase shifts which are quotients of Gamma functions, 
like (1.2). First note that in the case of the weight w given by (3.3), one has, from 



(3.22), (3.24), and (3.4), (3.5) 



Ki = bq, 9i = b,K2 = aq, 62 = a 
Ci = d/q,C2 = b/q,03 = 1/6, C4 = 1/a 



In order to recover ( 1.2 ) as a limit of the phase shift associated with the weight 
w we make the substitution >1 — e, a— 7>l + ea, 6— >'l + e6 where a and 6 are 
real numbers, before letting e — > 0. The final system is then 
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t{l-t)-- = tv(u + b){u + a)-t{t-l)v~^{u~a+l){u-h+l) 



dt 

dv 

t(l-t)— = v[2u + 2t{a-b-l)+a-b + 3] + 

2u{t -I) + a + b + 2 + t{a + b) + 
t;"^ [{2u - a)t + 2 + 3a + 6] 
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